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Abstract — In this paper a method of spectrum estimation of 
group delay using Auto- Regressive (AR) model is proposed. 
The method described is based on the properties of negative 
derivative of Fourier Transform (FT) phase function. The 
results obtained from the proposed method called modified 
product group delay are compared with the, Traditional group 
delay proposed by Hema.A Murthy proposedin 1992, modified 
group delay proposed by Rajesh M.Hegde in 2007 and product 
group delay proposed by Donglai Zhu and Kuldip K. Paliwal in 
2004 . This proposed method suppresses the spikes generated 
due to noise in spectrum and also provides better resolution. In 
the estimates obtained using this method the resolution 
properties of the periodogram estimate are preserved while the 
variance is reduced. This method works for even high noise 
levels (SNR = 0 dB or less). 

Index Terms — Auto-Regressive method (AR Method), 
group delay function, Product group delay, Spectral 
Estimation. 


I. INTRODUCTION 

The objective of this chapter is to explore an approach to 
spectrum estimation from the Fourier transform phase of the 
signal. The method described is based on the properties of the 
second order derivative of Fourier Transform phase function. 
Various attempts [l]-[2] have been made to demonstrate the 
spectrum estimation based on the properties of the negative 
derivative of the FT phase function, also called group delay 
function. The most important properties of the group delay 
function are the additive and high resolution properties [3]. 
Here the resolution refers to the sharpness of the peaks in 
group delay function, which is due to FT magnitude function 
behavior of the group delay function. 

Traditionally, the phase spectrum of the signal has been 
ignored, primarily because only the principal values of the 
phase can be estimated from the Fourier transform. For the 
phase to be used, the phase function will have to be 
unwrapped to produce a continuous estimate [3]. On the other 
hand, the group delay function [6] (defined as the negative 
derivative of the phase function), which has properties similar 
to the phase, can be computed directly from the signal. 

The most important properties of the group delay function 
are the additive and high resolution properties. Here 
resolution refers to the sharpness of the peaks in the group 
delay function which is due to the squared Fourier transform 
magnitude function behavior of the group delay function near 
the peaks. The key idea behind the new spectrum estimation 
method is that the properties of the group delay functions for 
noise and an auto regressive process are distinct. Group delay 
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is an important feature of the signal that can help in enhancing 
the signal quality in noisy conditions. Previous research 
works have revealed the usefulness of group delay in many 
applications. Recent studies on speech perception have 
revealed the importance of the phase of speech Signal. In 
order to overcome the problem of spikes in group delay, some 
researchers have suggested solutions such as modified group 
delay and product spectrum[4]-[5]. 


II. Group Delay 


The Group delay is defined as the negative derivative of the 
Fourier transform of a signal. Computing of group delay is 
difficult due to wrapping of the phase function [7]. 
Calculation of Group Delay is being done as follows. 

a) 

The X(bJ ) as a function of phase and magnitude can be 
expressed as, 

X{oi) = \X{0i)\e jg> ^ 

Now the group delay function can be expressed as 

f \ -d0(<u) 

t{oj) = 


UQJ 


( _ X R (t3)Y R ((D)+X I ( ( 3)Y](ts) 


( 2 ) 


where X R is Real part of the Fourier Transform of x(n), X j is 
Imaginary part of the Fourier Transform of x(n), Y R is Real 


part of the Fourier Transform of nx(n) and Y { is Imaginary 

part of the Fourier Transform of nx(n) 

2.1 Autoregressive Process in Noise (Estimation of the AR 
Spectrum) 

JlTj (n) = s(n) + u(n ) 


4 

s(n ) = - Y a k s(n -k) + Ge(n) (3) 

k= 1 

Where the excitation e(n) is a white Gaussian noise of 
variance unity and u(n) is an additive noise with variance 
dependent upon the coefficients are:ai=-2.760, a 2 = 3.809, a3 
= -2654, and a 4 =0,924. And the signal used is a two sinusoids 
in noise (estimation of frequencies of the sinusoids) 

x 2 in) = VlO exp [ j27i(0. 10 )n\ + V20 exp [j'2;r(0. 15)n\ + u(n) (4) 


Where u(n) is additive white Gaussian noise with the variance 
dependent upon the SNR. These examples are similar to the 
ones used in [8] for discussion of periodogram estimates. 

We assume a sampling frequency of 10 kHz and number 
of samples N=512 for AR spectrum estimation(4). Different 
realizations of x^n) and x 2 (n) are obtained by using different 
noise sequence each time. Spectral estimation calculated and 
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presented by using (3) , different method like MGDF, and 
product group delay function are compared here. 

III. Product Group Delay: 


Where S(&>) is cepstrally smoothed spectrum of X(ft>) in 
order to reduce spikes in group delay function. The parameter 
A and B are varied accordingly and estimated the spectrum . 


In order to make the GDF meaningful, a modification to the 
GDF has been proposed by replacing the power spectrum 

i i 2 

with the cepstrally smoothed power spectrum 

\S(<o)\ in (1). This gives the Modified Group Delay 
Function (MGDF) as follows [4]. 


<«> = 


__ X R (a)Y; (» 


(5) 


ls(«)l a 

We define the product spectrum Q(co) as the product of the 
power spectrum and the GDF as follows 


Q{co) = X(a>) 



= X r (oj)Y r (co) + X i (oj)Y i (co) 



The product spectrum [5] is influenced by both the 
magnitude spectrum and the phase spectrum. It enhances the 
region at the formants over the Modified Group Delay 
Function and has an envelope comparable to that of the power 
spectrum. 

In order to prevent the spikes on the group delay of signal, 
we will use a modified group delay function with different 
parameters as follows: 


T 



X R ( fij)Y fl tflj) TJfj l> )Y] l» 

S(Vr r ) " a 



( CO ) = t(co)\t(co) 


b -1 


(7) 

( 8 ) 


IV. Results: 

The results of the group delay function, modified group 
delay function product group delay and modified group delay 
function are presented for different values of the 
parameters a and b .The Fig. 1(a), 1(b) and 1(c) are 50 
overlaid , single realization and averaged spectrum estimation 
of AR Process respectively for the parameter a=0.9.similar 
results are depicted in Fig. 1(d), 1(e) and l(f).From these , 
average spectrum is giving a smoothed spectrum estimation. 

The Fig. 2(a), 2(b) and 2(c) are 50 overlaid , single 
realization and averaged spectrum estimation of AR Process 
for modified group delay method respectively for the 
parameter a=l and b=l .Similar results are illustrated in Fig. 
2(d), 2(e) and 2(f). 

The spectrum estimation of Modified group delay with 
a=0.9 ,b=0.4 and a=l,b=0.4 has been depicted in the Fig. 3(a), 
3(b) and 3(c) are 50 overlaid , single realization and 
averaged spectrum estimation of AR Process respectively for 
the parameter a=0.9. similar results are depicted in 
Fig.3(d),3(e) and 3(f).Based on the estimations carried out it 
is observed that the modified group delay method is giving 
better and smooth spectral characteristics when compared 
with the traditional group delay[10] and product group delay 
method. 
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Fig.l Spectrum estimation product group delay of AR process for A=0.9 and 0.4 respectively 
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Fig2 spectrum estimation modified product group delay of AR process for A=1 and b=l respectively 
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Fig 3 spectrum estimation modified product group delay of AR process for a=9.0 ,b=0.4 and a=l, b=0.4 

Respectively 


V. Conclusion: 
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